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Abstract 

If the constituent particles of fluid phases interact via long-ranged van der Waals forces, 
the effective Hamiltonian for interfaces between such fluid phases contains - in lateral Fourier 
space - non-analytic terms ~ q 4 lnq. Similar non-analytic terms characterize the effective 
Hamiltonian for two interacting interfaces which can emerge between the three possible co- 
existing fluid phases in binary liquid mixtures. This is in contrast with the structure of 
the phenomenological Helfrich Hamiltonian for membranes which does not contain such non- 
analytic terms. We show that under favorable conditions for the bulk densities characterizing 
a binary liquid mixture and for the long-ranged interparticle interactions the corresponding 
effective Hamiltonian for a model fluid membrane does not exhibit such non-analytic contri- 
butions. We discuss the properties of the resulting effective Hamiltonian, with a particular 
emphasis on the influence of the long range of the interactions on the coefficient of the bending 
rigidity. 



1 



1 Introduction 



In order to be able to describe nonplanar configurations of interfaces and membranes, the 
derivation and use of corresponding effective Hamiltonians has been studied intensively [1-3]. 
Depending on the environment and their internal composition interfaces and membranes can 
display rather complex behaviors [4]. A particular class of such systems is formed by the 
ubiquitous fluid- fluid interfaces and fluid membranes. In the case of interfaces the effective 
Hamiltonian takes on a capillary-wave like structure [5] while membranes are usually described 
in terms of the so-called Helfrich Hamiltonian [6] . 

On the phenomenological level the effective Hamiltonian contains two types of contributions: 
the first is related to the possible change of the interface or membrane area and is controlled 
by the coefficient ao of the surface tension while the second contribution is proportional to 
the square of the local mean curvature of the interface or membrane and is controlled by the 
coefficient k of the bending rigidity. In the following we consider fluctuating interfaces or 
membranes which are planar on the average and do not change their topology; thus contribu- 
tions due to the Gaussian curvature do not matter. In lateral Fourier space the contribution 
from the g-mode /(q) of a local height configuration to the effective Hamiltonian is propor- 
tional to |/(q)| 2 (qa) 2 cr(q), where a is a microscopic length scale proportional to the particle 
diameter and a(q — > 0) = ao + k (qa) 2 . 

Here we focus on the ensuing structure of the effective Hamiltonian for systems in which the 
interparticle interactions are of the long-ranged van der Waals type. This issue becomes acute 
if one tries to justify and to derive the phenomenological capillary-wave Hamiltonian from a 
microscopic theory such as, e.g., density functional theory [5]. In such approaches it turns 
out that for interfaces between fluid phases in systems governed by long-ranged forces the 
effective surface tension a(q) exhibits the form a(q — > 0) = do + o\ (qa) 2 ln(qa) + n(qa) 2 , 
and thus contains a leading non-analytic term a± (qa) 2 ln(qa) with <j\ > which is not cap- 
tured by phenomenological approaches. This logarithmic singularity in Fourier space can be 
traced back to the divergence of the third and higher moments of the interparticle interac- 
tion potentials decaying as function of the distance ~ r -6 . For fluid interfaces the presence 
of such a non-analytic contribution has been established theoretically [7-10]. This implies 
that for small g ^ one has a(q) < ao which has been confirmed also experimentally for 
various systems [11] as well as in some simulations [12] but not in all [13]. On the other 
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hand such non-analytic terms are absent in the effective Helfrich Hamiltonian for membranes 
which, however, successfully describes various properties of fluid membranes. This is puzzling 
because the particles making up membranes invariably also exhibit long-ranged van der Waals 
interactions which in turn should lead to non-analytical bending contributions. 
Our objective is to construct a simple model of a fluid membrane based on the extension of 
a model of two interacting fluid-fluid interfaces. We want to check under which conditions, if 
any, the absence of non-analytic terms of the type 0\ (qa) 2 ln(qa) in the effective Hamiltonian 
for a membrane is possible, and what kind of influence on the remaining terms these condi- 
tions have. In the following section we recall the relevant results concerning the structure of 
the capillary- wave Hamiltonian. In Sect. [3] we discuss a simple model of fluid membranes 
in a system with long-ranged forces which is based on a model of two interacting fluid-fluid 
interfaces in a binary liquid mixture. We establish the conditions under which the effective 
Hamiltonian for the fluid membrane is free from non-analyticities present in the correspond- 
ing capillary-wave Hamiltonian for the interface and is thus compatible with the structure of 
the Helfrich Hamiltonian. In Sect. 0] we compare our predictions for the resulting effective 
Hamiltonian with those discussed in the literature. 

2 Effective Hamiltonian for a fluid-fluid interface 

In this section we recall the basic facts pertinent to the structure of the capillary-wave Hamil- 
tonian H C w[f] for a fluid-fluid interface. Its local height relative to the reference plane 
z = is described by the function z = /(R), where R = (x,y) denotes the lateral coor- 
dinates. Various aspects of this structure have been discussed in the literature. In particular, 
the issue of a local versus a non-local structure of Ttcw[f] has been extensively analyzed 
for the cases of short-ranged (exponentially) and long-ranged (algebraicly) decaying interac- 
tions [7-10, 14]. A suitable framework for analyzing such issues is density functional the- 
ory for non- uniform fluids. This analysis is particularly straightforward if the non- uniform 
one-component fluid density p(R, z) associated with an interface configuration /(R) is ap- 
proximated within the so-called sharp-kink approximation by a piecewise constant function 
p s hfc(R, z) = pp 0(/(R) — z) + p a @(z — /(R)), where p a and pp denote the bulk densities of 
the coexisting fluid phases a and /?, and 0(z) denotes the Heaviside function. If, moreover, 
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the effective Hamiltonian is truncated to be bilinear in /, it can be written as [7-9] 

n cw [f] = \ J ^ |/(q)|%M<?) (i) 

where 

/(q) = y^/(R)exp(-iq-R) . (2) 

The wavevector dependent surface tension a(q) in Eq.JJJ is given by 

o-(q) = q~ 2 [w(q) - w(0) } (p a - ppf , (3) 

where w{q) denotes the Fourier transform of the long-ranged part of the spherically symmetric 
interparticle interaction potential w{r = |(R, z)\) taken with respect to the lateral coordinates 
for z = 0: 

w(q) = J d 2 Rw(\(R,z = 0)|) exp(-iq-R) . (4) 

It has turned out to be suitable to adopt for the long-ranged part of the van der Waals pair 
potential w(R, z) the form 

m(R^)i) = - (R2+z 1 +a2)3 , (5) 

where a corresponds to the hard core radius of the fluid particles and A > characterizes the 
strength of the attractive interparticle interaction. For q = qa 1 the ensuing a(q) has the 
following non-analytic form: 

a(q) = a + a x f ln{q) +a 2 <f + O(cf) , (6) 

where a = £g (p a - p p ) 2 > 0, a x = \a Q > 0, a 2 = \ a C < 0, C = C s - 3/4 - In 2 = 
—0.866, and Cs denotes Euler's constant. 

A more realistic approach to determine a(q) [9] takes into account the influence of local 
interfacial curvatures on the actual smooth intrinsic density profile. The effective Hamiltonians 
for interfaces both in one-component [9] and in binary liquid mixtures [10] have been analyzed 
along these lines. For long-ranged van der Waals interactions in each case the presence of non- 
analytic terms in a(q) (Eq.JT])) has been established. 
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3 A model of a fluid membrane 



For the comparison between effective Hamiltonians for fluid-fluid interfaces and fluid mem- 
branes it is particularly suitable to consider binary liquid mixtures. Upon special choices of 
the thermodynamic conditions these systems allow for the coexistence of three fluid phases 
denoted as a, (3, and 7. In the presence of appropriately chosen boundary conditions or ex- 
ternal ordering fields one can consider a situation in which a layer of - say - phase f3 with 
mean thickness £ separates the phases a and 7 [15]. In such a system there are two fluid-fluid 
interfaces the positions of which are denoted by / Q( g(R) and I + //3 7 (R). They separate the 
phases a, (3 and /3,7, respectively, (see Fig. 1). We note that although the system is char- 



Figure 1: The system under consideration consists of two fluid-fluid interfaces f a/ 3(H) and ^+/g 7 (R) 
separating the phases a, (3 and (3, 7, respectively, and fluctuating around their mean positions z = 
and z — £, respectively. 

acterized by six number densities pi\, with i = 1,2, and A = a,f3, 7, where pi\ denotes the 
number density of the i-th component in phase A, three-phase coexistence allows for only one 
independent thermodynamic variable such as temperature; on the corresponding triple line 
the chemical potentials pi(T) and P2(T) of the two species are fixed. In addition there are 
three interparticle interactions present in the system: two among the two species and one 
between the different species. They are assumed to be spherically symmetric and are denoted 
by Wij(r) = Wji(r) = Wy(|(R,z)|) with i,j = 1,2. 

For such a system containing two interfaces the capillary-wave Hamiltonian 
T^cto ( [/a/3 ) /ft"/], £) is a functional of the two interfacial positions / Qj a(R) and //3 7 (R) and a 
function of the distance £. Applying the sharp-kink approximation to the density functional 
for binary liquid mixtures described, e.g., in Refs. [8,10] yields within the bilinear approxima- 
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tion, which retains non-locality, the following form 

1 f d 2 q 



+ 



where 



1 f Uj 

HcwitfapJihht) = 2 J (2^2 { 2 /a/3(q)//3 7 (-q) w af3 ^(q,£) (7) 

2 |2 
/a/?(q) [^^(g) - *q/3,/3 7 (0,^)] + //3 7 (q) [^ 2 0-/3 7 (9) ~ *a/3,/3 7 ( ^) 



2 /• 

w\ K>r) s{q, z) = ^ (pa - Pin)(Pjr, ~ Pjs) / d 2 -R exp(-iq-R)u7jj(R,2r) 



(8) 



and [8,10] 



0"A«(<?) = Q 2 [™A K ,A K (<?, 2 = 0)- w\ Ky \ K (q = 0,z = 0)] 



(9) 



The above results can serve as a starting point to construct a simple model of a fluid membrane. 
To this end we take the two interface configurations to be in phase, i.e., / aj a(R) = //3 7 (R)- 
This renders a model fluid membrane consisting of phase (3 embedded on one side by phase a 
and on the other side by phase 7. The thickness t of the membrane is uniform and its upper 
and lower boundaries have the same shape described by /(R) = / Q/ g(R) = /^(R). In this 
case and within the bilinear approximation the capillary-wave Hamiltonian reduces to: 
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«([/],<) = - 2 j (2 , )2 



/(q) 7(<7)<7 2 , 



(10) 



with 



Q 2 1(q) = {[ (Pia - Pip)(Pja ~ Pjp) + (Pi/3 - Pij)(pjp - Pj-y)] [u>ij(q, 0) - 7X^(0, 0)] 

+ 2(p ia - Pip){pjp - p h )[wij{q,t) - Wij{0,£)] } . (11) 



With the choice 



Wij (\(R,z)\) = - 



(R2 + z 2 + a 2. )3 

for the long-ranged interparticle potentials one has 



(12) 



Wij(q,£) - Wij(0,£) 



+ ^ 2 



+ ^ (ln(9>/4 + + C °) 



(13) 



For reasons of simplicity in the following we assume Ojj = a, i,j = 1,2. This choice leads to 
the following expression for 7 (q): 



7(9) = 70 + 71 q 2 In 9 + 72 <7 2 + CKq 4 ) , 



(14) 



where, with £ = £/a, 
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+ - — - 2 (Pi« - Pip){PjP ~ Pit) | > ( 15 ) 

2 

71 = X] A V (P™ ~ Pij)(Pja ~ Pj-y) , (16) 

and 

2 

72W = ^2 A v { ln ( J + ^) (P ia ~ Pip)(PiP ~ Ph) 

+ ^0 (pi Q - Pi-,){pj a ~ Pj-y)} ■ (17) 

As expected, similar to the case of single interface (Eq. (j6]l ) the effective Hamiltonian for 
the model fluid membrane contains a non-analytic contribution 71 q 4 ln q. In this sense the 
structure of the effective Hamiltonian given by Eqs, (|10|14flT7]) is not compatible with the 
phenomenological Helfrich Hamiltonian ansatz which for small membrane ondulations can be 
expressed in its form as in Eq. fTT4|) but with 71 = 0. 



Our purpose is thus to find conditions under which the coefficient 71 of the non-analytic 
contribution in Eq. (ll4p vanishes. There are two particularly simple choices of the number 
densities pi\ and the amplitudes A^ of the interaction potentials which fulfill this requirement. 
The first choice (I) puts constraints on the densities of the phases a and 7 and stipulates 



Pi 



Pi-y 



i=l,2 (I). 



This condition imposes that the two phases on both sides of the membrane are identical. The 
second choice (II) puts constraints both on the interaction amplitudes and on the densities. 
First, it requires that 



A 12 = V A n A 22 [He 



(19) 



which leads to 



71 



7T 

32 a 2 



^4.11 (Pla - Pl-y) + \A422 (P2a ~ P2 7 ) 



(20) 
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The additional requirement 



An (Pia ~ Pl-y) 



A 2 2 (P2a ~ P2-y) (I lb) 



(21) 



implies 71 = 0. It is straightforward to show that the above condition I (Eq. (fl~8l) ) leads to 



7T 

4^2 



r An(p Xa - p lf3 ) + \[A~2~ 2 (p 2(x - P2/3) 



(22) 



where 



1 + 



(23) 



Interestingly, if the conditions (Ha) (Eq.JT9|)) and (lib) (Eq. ([2T1) ) are imposed, the corre- 
sponding effective surface tension 7^(9) has exactly the same form as for the first condition, 
i.e., j( H \q) = ^^(q)- The fact that the requirements (//), which put constraints on both 
the densities and the interaction amplitudes, lead to the same result as the requirement (I), 
which identifies the phases a and 7 but does not involve the interaction amplitudes Aij, can be 
understood as follows. We consider a typical contribution to the free-energy density functional 
which describes the interaction between particles located in a region V Q of the binary liquid 
mixture with a specific particle of type k, k = 1,2, located at r' somewhere in the system. 
This term is proportional to 



/ 'Pr V] w ik (r - r')p ia ~ V" A ik p ia = A lk p la + A 2k P2c 

JVr, : 



(24) 



'A 



kk- 



An Pia + \TAnp2a) = V Akk ( V^4iiPi7 + \fM2P2-A = A lk pi y + A 2k p2j 



where the conditions in Eq. (fT9|) and Eq. (|2T|) have been used. One concludes that this con- 
tribution to the free-energy functional has the same form as if the region V a would be filled 
with particles with densities pi-y instead of pi a . But this is exactly the requirement in Eq. (fT8l) 
corresponding to choice (I) which identifies the phases a and 7. 
In the next section we discuss the properties of the resulting effective Hamiltonian. 



4 Discussion 

In the previous section we showed that for special choices for the densities or the interparticle 
interactions in binary liquid mixtures there are no non-analytic contributions to j(q) in the 
limit of small q (up to and including 0(q 2 ). This choice eliminates the leading non-analytic 
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contribution for any membrane thickness £, because 71 does not depend on £ (see Eq.(fl6j)). 
It turns out that independent of whether constraints of type (I) in Eq. (fT8j) or of type (II) 
in Eqs, (fT9l I2T1) are imposed the resulting effective Hamiltonian for the model fluid membrane 
takes the form given by Eqs.(HQD and {22]). The function 7 (/) (<?) in is determined 

by the bulk number densities pi\,i = 1,2, X = a, (3, the interaction strengths A\\,A22, the 
particles size a, and the membrane thickness £. The function 7^(9) factorizes into a product 
of two functions. The first factor depends on the densities and interaction strengths only and 
is non-negative. The second factor depends on q and parametrically on £ only; the parameter 
a sets the scale for the variables q and i. This second factor, which we denoted as <p(q,£), is 
particularly interesting because - contrary to the first factor - it can change sign depending on 
the values of q and £. This possibility of 7^ (q) to change sign appears because the coefficient 
72 in Eq. lfTil) is inherently negative, i.e., the contribution from the long-ranged forces to the 
coefficient of the bending rigidity is negative. (Note that within the sharp-kink approximation, 
which takes only the influence of the long-ranged forces into account, also 02 is negative (see 
Eq.([6]) and the expressions following it).) This conclusion checks qualitatively with a recent 
analysis by Dean and Horgan [16] who have expressed the coefficient of the bending rigidity 
in terms of the membrane thickness and the dielectric constants of the membrane (e) and of 
the surrounding medium (e'): 



{l) = -vk(—) n i+r J • < 26 > 

Dean and Horgan [16] have not discussed the issue of the presence of non-analytic terms in the 
effective Hamiltonian. However, our result and those in [16] agree concerning the functional 
form of the dependence of the coefficient 72 (£) on the membrane thickness £. Within both 
approaches the coefficient of the bending rigidity depends logarithmically on the membrane 
thickness, i.e., 72 ~ ln(l + £ ). Of course realistic membrane models yield additional contribu- 
tions to the bending rigidity stemming from other types of interactions present in the system. 
In our approach only the long-ranged contributions to the bending rigidity are considered. 
In this latter case the negative coefficient 72 of the bending rigidity in the presence of the 
positive coefficient 70 of the surface tension leads to an instability at small wavelengths of the 
membrane ondulations. According to Eq. ([23l) this instability occurs for 

8 £ 2 

q 2 > =T5 =75- (unstable; I, II) . (26) 

(1 + f) ln(l + f ) 



On the other hand the wavevectors must be smaller than the physically allowed maximal one 
Qmax i$ 1- This implies that the values of i for which the instability can occur fulfill the 
condition 

£ > E = exp (-^-) (unstable; I, II) . (27) 

V Qrnax J 

For q m ax = 1/2 one has £o = 9 x 10 6 . This condition states that for membrane thicknesses 
< £ < £q the negative bending rigidity coefficient does not give rise to instabilities for mem- 
brane ondulations with wavevectors within the physically accessible range q < q m ax- 



Finally we mention that the vanishing of the coefficient 71 can also occur in binary liquid 
mixtures in which the interactions wn and W22 are repulsive, i.e., An, A22 < while the 
interactions w±2 are attractive, i.e., A12 > 0. (It is conceivable that such a situation may arise 
in multicomponent complex fluids with effective interactions between two dominating species 
upon integrating out the degrees of freedom of the smaller species. This can occur if the two 
species are oppositely charged.) This is a different situation from the one considered above in 
which all long-ranged interactions were assumed to be attractive, i.e., An , A22 , A\2 > 0. In 
this present case the conditions (Ha) and (lib) are replaced by 



A 12 = V(-Aii)(-A 22 ) (IHa) (28) 



so that 



7T 

71 = "32^2 



V~ A U (Pla - Pl 7 ) - V~ A 22 (plot ~ P2 7 ) , (29) 



and by 



y/-Au ( Pla - pi 7 ) = (P2a ~ P2 7 ) (J I lb) , (30) 



respectively. It is straightforward to see that in this case the effective surface tension denoted 
as ^ In ^ (q) is given by 

J (III) (q) = -£f [ V^MPIc* - Pip) - V ZA 2~2(p2a - P2/3)] <p{q,?) ■ (31) 

Accordingly the model fluid membrane is unstable with respect to long-wavelength ondula- 
tions: 

8 £ 2 

q 2 < =n =2~ (unstable; III) . (32) 



(1 + £ 2 ) ln(l + l 2x 
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This implies that for ondulations with given g-values only membranes with thicknesses 



are unstable. 

To summarize, we have shown that it is possible to choose conditions under which the leading 
non-analytic contribution to the effective Hamiltonian of a fluid membrane in the presence of 
long-ranged forces vanishes. One of them amounts to the requirement that the embedding 
phases on both sides of a fluid membrane are identical. We have also checked that the con- 
tribution from long-ranged forces to the coefficient of the bending rigidity is negative and we 
have discussed the implications on the stability of membranes with respect to ondulations. 
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